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Abstract 

We present the result for the invariant Bk factor of K°-K° mixing in the chiral limit and to 
next-to- leading order in the 1/N C expansion. We explicitly demonstrate the cancellation of the 
renormalization scale and scheme dependences between short- and long-distance contributions in 
the final expression. Numerical estimates are then given, by taking into account increasingly re- 
fined short- and long-distance constraints of the underlying QCD Green's function which governs 
the Bk factor. 



lr This is a revised version of our article in Phys. Lett. B490 (2000) 213. 



1 Introduction 



The Standard Model predicts strangeness changing transitions with AS = 2 via two virtual W- 
exchanges between quark lines, the so-called box diagrams. The low-energy physics of these transi- 
tions is governed by an effective Hamiltonian which is proportional to the local four-quark operator, 
(summation over colour indices within brackets is understood and <?l = 



q, 



!AS- 



-_ 2 (x) = (s L (x)j' 1 d L (x))(s L (x)j IJ ,d L (x)), 



(1.1) 



modulated by a quadratic form of the flavour mixing matrix elements X q = V* d V qs , q = u, c, t , with 
coefficient functions i*i,2,3 of the heavy masses of the fields t, Z°, W^, b, and c which have been 
integrated out 0: 



n,AS=2 



G 2 F M^ 
4tt 2 



[A*Fi + X 2 t F 2 + 2X c X t F 3 ] Cas=2(h)Qas=2(x) ■ 



(1.2) 



The operator Qas=2 is multiplicatively renormalizable and has an anomalous dimension "f(a s ) which 
in perturbative QCD (pQCD) is defined by the equation 



~j n < Qas=2 >= --7 «s < Qas=2 > , l(a s ) = — 7i + — 



72 



with (/?! = ±(-11^ + 2^)) 

1 



7i 



1 - 



and 72 = 



32 



1 - 



19 



(1.3) 



(1.4) 



where k = in the naive dimensional renormalization scheme (NDR) and k — —4 in the 't Hooft- 
Veltman renormalization scheme (HV). The renormalization /i-scale dependence of the Wilson coef- 
ficient Cas=2(m) in Eq. ( |l]| ) is then, (/3 2 = ~j^N 2 + \^k^ n f + T5 N c n f) > 



Cas=2(m) 



72 



a B (^) /1433 



7i 
-A 



«s(M) 



-ft 



1 



1936 8 



JL i 

11 n c 



(1.5) 
(1.6) 



where the second line gives the result to next -to-leading order in the 1/N C expansion, which is the 
approximation at which we shall be working here. 
The matrix element 

< ^°|QAs= 2 (0)|if° >= p 2 K M 2 K B K (v), (1.7) 

defines the so-called Bk~ parameter of K° — K° mixing at short-distances, which is one of the crucial 
parameters in the phenomenological studies of CP-violation in the Standard Model. In the large--/V c 
limit of QCD, the four quark operator Qas=2 factorizes into a product of two current operators. Each 
of these currents, to lowest order in chiral perturbation theory, has a simple bosonic realization: 

(sl^l) = =► ^F hr [X Z2 {D^)U] = *F ff>K° + ■■■, (1.8) 

dln{x)z 2 2 V2 

where Fq, denotes the coupling constant of the pion in the chiral limit and U is the 3x3 unitary matrix 
in flavour space, U = 1 + • • •, which collects the Goldstone fields and which under chiral rotations 
transforms like U — > VrUV^. The large-iV e approximation, with inclusion of the chiral corrections in 
the factorized contribution, leads to the result: 

3 

B k \n c ^oo = j and Cas=2(h)\n c ->oo = 1 • (1-9) 

2 For a detailed discussion, see e.g. refs. ly, H and references therein. 
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In full generality, the bosonization of the four-quark operator Qas=2(x) to lowest order in the chiral 
expansion is described by an effective operator which is of 0(p 2 ) ^: 

F 4 

QasMx) => — ^ 3as= 2 (m) tr [X 32 (D^)U\ 32 (D^)U] , (1.10) 

where A32 denotes the matrix A32 = Si 3 8 2 j in flavour space and <7as=2(m) is a dimensionless coupling 
constant which depends on the underlying dynamics of spontaneous chiral symmetry breaking (S%SB) 
in Q CD . The relation between the coupling <?AS=2(/i) and the phenomenological Bk factor defined in 
Eq. ([Tt]) is simply 

3 

B K (ii,m Utd ,s -* 0) = -5as=2(m) • 
Notice that the coupling g&s=2{n), and hence Bk, is perfectly well defined in the chiral limit, while 



the matrix element in Eq. (1.7) vanishes in the chiral limit as a chiral power. To lowest order in the 



chiral expansion, the relation to the so called invariant £?A'~factor is as follows: 

3 

Bk = ^Cas= 2 (^) x 9 , as=2(a*) , (1.12) 

which means that the coupling gAS=2 (m) must have a /1 dependence (and a scheme dependence) which 
must cancel with the \x dependence (and scheme dependence) of the Wilson coefficient Cas=2{^)- 
The purpose of this note is to show how the mechanism of cancellations works in practice within the 
framework of the 1/N C expansion. This will allow us to give a numerical result for Bk which is valid 
to lowest order in the chiral expansion and to next-to-leading order in the 1 /N c expansion. We wish 
to emphasize that the novel feature of this work is that, to our knowledge, it is the first calculation of 
Bk which explicitly shows the cancellations of scale and scheme dependences. 



2 Bosonization of Four— Quark Operators 

The QCD Lagrangian in the presence of external chiral sources l^, of left- and right- currents, but 
with neglect of scalar and pseudoscalar sources which is justified in the chiral limit approximation we 
shall be working here, has the form 



C 



QCD 



(2.1) 



with T> x the Dirac operator 



2>, 



1-75 



1 



'75 



2 ^2 

The bosonization of the four-quark operator Qas=2{x) is formally defined by the functional integral 

(QAS=2(x)) - 

- / *v 



(2.2) 



Tr V-H-^c 
x *«„(*) 



32 



d Q -iq-jx-y) 
(27T) 4 



Tr V-H 5Vx 



T) _1 7 ±. 

x 8l»(y) 



:',2 



(2.3) 



where the trace Tr here also includes the functional integration over gluons in large N c . The first 
term corresponds to the factorized pattern and gives the contributions of 0(N 2 ); the second term 
corresponds to the unfactorized pattern and it involves an integral, (which is regularization dependent,) 
over all virtual momenta q. This is the term which gives the next-to-leading 0(N C ) contribution we 
are interested in. 



3 See e.g. refs. 
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To proceed further, it is convenient to use the Schwinger's operator formalism. With /P the 
conjugate momentum operator, in the absence of external chiral fields, the full quark propagator is 



1 i 

(x\—\y) = (x\ ^ rx — rly) with {Af\y)=i tJ 



i-^. Cm 

On 



5{x - y) . (2.4) 



The chiral expansion is then defined as an expansion in the l a and r a external sources. The precise 
relation between the formal bosonization in Eq. (2.3) and the explicit chiral realization in Eq. ( l.lCj ) 
can best be seen from the fact that: D^U^U — d^Wlf + iWr^U — il^. This shows that there are 
a priori six different ways to compute the constant QAS^ifJ-), although they are all related by chiral 
gauge invariance. One possible choice is the term 



4z =2 (z) = ~ 5A5= 2 (^) {• ■ ■ -tr[A32t/Vt/A 3 2£/ + 7>[/] + ■ • •} . (2.5) 

This is a convenient choice because the underlying QCD Green's function is the four-point function, 
[ ^ d (x) ee E^fr"^^) and R^ix) ee J2j^ x )^l^ s a (x),] 

W££&(?.0 = I™ ^ ^ I d^x dS d A z e*' e~^ (0\T{L» d (x) R%{y) L^(0) <(z)}|0) , (2.6) 



In fact, what we need, as seen in Eq. (2.3), is the integral of the unfactorized four-point function 
Wlrlr( ( 1-> over t ne four-vector q with the Lorentz indices of the two left-currents contracted. This 
is a quantity which is a good order parameter of SxSB. The integral over the solid angle has the form, 
(Q 2 = ~q 2 ,) 

dn q gu.WrnUq, 01 unfactorized = " 9° P ) ^LRL R (Q 2 ) . W 

where the transversality in the four-vector / follows from current algebra Ward identities. We are still 
left with an integral of the invariant function VVj j 1 ^ LH (Q 2 ) over the full euclidean range: < Q 2 < oo 
which has to be done in the same renormalization scheme as the calculation of the short-distance 
Wilson coefficient Cas=2( x ) m Eq. (1.2) has been done, i.e. in the MS-scheme. The coupling 
constant <?as=2(a0 defined with dimensional regularization in d = 4 — e dimensions is then given by 
the following integral, 

where, for convenience, we have normalized Q 2 to a characteristic hadronic scale ^ nad , (z ee Q 2 //i nad ,) 
which in practice we shall choose within the range: 1 GeV < //had. < 2 GeV. 



3 Bk to Next— to— Leading Order in the 1/N C Expansion 

In full generality, Green's functions in QCD at large N c j^] are given by sums over an infinite number 
of hadronic poles M. Since the function W^ ] RLR (Q 2 ) is an order parameter of spontaneous chiral 
symmetry breaking it receives no contribution from the perturbative continuum and satisfies an un- 
subtracted dispersion relation. After repeated use of partial fractions decomposition, one can see that 



for the particular case of the function W[z], [see Eqs. (2.8), (2.7) and fl2.q)] the most general ansatz 



ie70. 



in large- N c QCD is an infinite sum of poles, double poles and triple poles n of the form 

oo oo oo oo c 



4 Other choices are indeed possible but, in general, the underlying QCD Green's functions have pieces which are not 
order parameters and spurious contributions which depend on the regularization. We have shown the equivalence among 
the various possible choices, but we postpone the detailed discussion which is rather technical to a longer publication. 

5 We thank J. Bijnens and J. Prades for pointing out to us the existence of triple poles, which were ignored in the 
previous version of this manuscript. 
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where pi = Mf / and Mi is the mass of the i-th narrow hadronic state. Triple poles result from 
the combination of two odd-parity couplings involving vector mesons and Goldstone particles. The 
first term on the r.h.s. is the contribution from the Goldstone poles to the W [z] function. This 
constant term is known from previous calculations, (see refs. |Q |8[ |{|) with which we agree. Constant 
terms, however, do not contribute to the integral in Eq. fl2.gp defined in dimensional regularization. In 
the absence of a solution of large-iV c QCD, the actual values of the mass ratios pi, and residues on , 
Pi and 7i remain unknown. As discussed below, there are, however, short-distance and long-distance 
constraints that the function W[z\ has to obey. With W[z] given by Eq. (3.1), the integral in Eq. (2 
can be done analytically, with the result 



0AS=2 (M) e ) = 1 



Mhad. 

32^ 2 F 2 



- + log 47T - 7jB + log /Lt 2 /Mhad. + 1 ) E ' 
6 ' i=l 

OO OO q ^ CO \ 

-£«iiogpi+53-i+-£^[ . 

i=l i=l PI z i=l p * J 



(3.2) 



We shall next explore the short-distance properties and long-distance properties which are pres- 
ently known about the function W[z], and hence about the parameters c^, /3j, ji and pi. 

• The large-Q 2 behaviour of the function W^ iH ((3 2 ) is governed by the OPE. We find the result 



7T 



1 



e 1 



Q 4 



(3.3) 



26 ( 5 + ^ + (v 

with k = in the NDR scheme and k = —4 in the HV scheme. The fact that the residue of the 
l 

Q 1 



L power term in the OPE is known, entails the constraint 



i=l 



OLi — 1Z + —S , 



with 

K = 



24ir 2 — + 

7T 



N r ai 



F 2 

Phad. 



and S = 



47r 2 (5 + K )^ + 



This allows us to define the renormalizcd coupling constant 5^5=2 (W ^ n ^ ne MS-scheme 



N n ai 



F 2 

Mhad. 



Mhad. 
327r 2 J F 2 



7llogM 2 /Mh a d.+^ + <5- 



E 

i=l 



<*s(m) 
a S (Mhad.) 



1 A^had. 

2 16tt 2 F 2 



ft + S + J] 



-ttj log p, 



-on log p 4 




(3.4) 



(3.5) 



(3.6) 



,(3.7) 



where in the second line we have written the renormalization group improved result. We insist 
on the fact that this result, contrary to the various large-7V c inspired calculations which have 
been published so far fu], y, |, % || § | is the full result to next-to-leading order in the 
1/N C expansion. The renormalization /i scale dependence as well as the scheme dependence in 
the factor S cancel exactly, at the next-to-leading log approximation, wit h the short distance 
p and scheme dependences in the Wilson coefficient Cas=2(p) m Eq. (1.6). Therefore, the full 
expression of the invariant -Br- factor to lowest order in the chiral expansion and to next-to- 
leading order in the 1/N C expansion which takes into account the hadronic contribution from 
light quarks below a mass scale /ihad. is then 



B 



K 



Ci s (p 



had. ; 



had. , 



1229 



1936 



O 



N c a 2 (p 



had. 



Mhad. 

32tt 2 F 2 



E 



-OLi l°g Pi H 1" 77-9 

Pi 2 pf 



(3.8) 
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The numerical choice of the hadronic scale //had. is, a priori, arbitrary. In practice, however, 
/ihad. has to be sufficiently large so as to make meaningful the truncated pQCD series in the 
first line. Our final error in the numerical evaluation of Bk will include the small effect of fixing 
A*had. within the range 1 GeV < /ihad. < 3 GeV. 

The fact that there is no 1/Q 2 term in the OPE expansion of the function VV| i 1 ^ Lfl (Q 2 ) implies 
the sum rule 

OO 00^,00 

i=l p% i=l p l i=l p l 

This sum rule plays an equivalent role to the 1st Weinberg sum rule in the case of the Hlr{Q 2 ) 
two-point function \ Like in this case, it guarantees the cancellation of quadratic divergences in 



the integral in Eq. (2.8) and it relates the contribution from Goldstone particles to a specific sum 
of resonance state contributions. What this means in practice is that, although the Goldstone 
pole does not contribute to the integral in Eq. ( |2.S| ) defined in dimensional regularization, the 
normalization of the function W[z] at z — is, precisely, fixed by the residue of the Goldstone 
pole contribution; i.e., W[0] = 6. 

The slope at the origin of the function W[z] is fixed by a linear combination of coupling constants 
of the 0(p 4 ) Gasser-Leutwyler Lagrangian with the result [] 

OO OO r. OO 2 

E^+ 2 E^+ 3 Et = 24%^ + 5L 2 + L 3 + L 9 ) . (3.10) 

tip* tip* t=itf F o 



The numerical evaluation of Bk in Eq. (3.8) requires the input of the mass ratios pi and couplings 



a.i, (3i and 7i of the narrow states which fill the details of the hadronic spectrum of light flavours. 



Our goal is to find the minimal hadronic ansatz of pole terms of the form shown in Eq. (3.1) which 



satisfies the short-distance and long-distance constraints discussed above, and use then this minimal 



hadronic ansatz to evaluate the integral in Eq. (2.8). The fact that the function W[z] is an order 



parameter of S^SB ensures that it has a smooth behaviour in the ultraviolet and it is then reasonable 
to approximate it by a few states. (There is no need here for an infinite number of narrow states to 
reproduce the asymptotic behaviour of pQCD, as it would be the case with a Green's function which 
is not an order parameter of SxSB.) This procedure has been successfully tested in other examples, 
like the calculation of the electroweak tt + — tt° mass difference |l3| and the decay of pseudoscalars 
into lepton pairs jL5j , and it has been applied |l6| to the evaluation of matrix elements of electroweak 
penguin operators as well Q 



Numerical Evaluation and Conclusions 



We now proceed to the numerical evaluation of the analytic result for Bk in Eq. ( p.8|) by successive 
improvements in the hadronic input. 

The crudest approximation is the one where the only terms of the next -to-leading order in the 
l/N c expansion which are retained are those from the first line in Eq. (3.8). This implicitly assumes 
that all the way down to the hadronic /ihad. scale there is only pQCD running and that the hadronic 
contribution below the /ihad. scale due to light quarks is negligible. One can see that from the integral 
in Eq. (|2.8|), if we write it in an equivalent cut-off form: 



dzW[z] + 



6 See e.g. ret@ 
7 The r.h.s. in Eq 
8 See also refs. 



r 2 ! fEjjP-L 



A 2 /V 



had. 



(4.1) 



agrees with the expression reported in ref. [^j . 



and refs 



0§ 



for other recent evaluations of electroweak penguin operators. 
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with A 2 = p? exp (l + ||), and where in the second integral we use the asymptotic OPE expression 
in Eq. ( |3.3| ). The approximation we are discussing is the one where the first low-energy integral is 
simply neglected. The result of this next-to-leading log (nil) approximation is to bring up the large- 
N c prediction in Eq. to a value which, including an estimate of higher order corrections, (we use 

(3) 

A^L = (372 ± 40) MeV,) but ignoring the systematic errors involved in the neglect of the hadronic 
contribution, would be 



B K \ 



0.96 ±0.04 



for 



AUiad 



= 1.4 GeV. 



(4.2) 



The problem, however, with this simple estimate is that the underlying assumption of neglecting 
entirely the low-energy hadronic integral does not satisfy any of the long-distance matching constraints 
discussed in the previous section. 

One can considerably improve on the previous estimate by taking into account the contribution 
from the lowest hadronic state to the low-energy hadronic integral, the p vector meson. It turns 
out that, in the presence of possible triple poles in Eq. ( |3.1| ), this is the minimal hadronic ansatz 
approximation at which one can fix the large- V c hadronic expression in Eq. (3.1) to satisfy the two 
matching constraints in Eqs. (pM and ([3.10]) as well as the first non-trivial constraint from the OPE 



in Eq. (3.3). In fact, it has been shown that the phenomenological values of the L,; constants which 
appear on the r.h.s. of Eq. (3.10) agree well with those obtained from the integration of the lowest 
vector state ^4], |2f| alone |j. In this approximation |^5) , the sum rule in Eq. ( p. 10 ) becomes simply 



a v py + 2(3 v p v + 3"f V 



21,4 , 



with pv — My I p 



had. 



(4.3) 



Using this constraint, together with Eqs. (3.4) and fl3.9| ) restricted to the lowest hadronic state, we 
can then solve for ay, Pv and 7v iu terms of pv with the result: 



a v = 2ATra s p v Fq/M v , (3 V = ~p v (3p v + 2a v ) , and j v = p v (9p v + a v ) . 



(4.4) 



The resulting W[z] function, normalized to its value at the origin W[0] = 6, is plotted in Fig. 1 below 
as a function of z = Q 2 / '/i 2 ad • Also shown in the same plot are the chiral perturbation behaviour of 
the W[z] function from the knowledge of its value and the slop e at the origin (the green line) and the 
corresponding behaviour from the OPE expression in Eq. ( |3.3[ ) (the blue line). We can now make an 
improved estimate of Bk by calculating the integral 



dzW[z] + 



a 2 /m 



had. 



dzW OPE [z} 



(4.5) 



with W[z] in the first integral approximated by the minimal hadronic ansatz just discussed. The choice 
of z which minimizes the dependence of Bk on z is the one at which the two curves W[z] and Wope[z] 
intersect. In terms of this z, which separates the long-distance part of the integral estimated with the 
minimal hadronic ansatz and the short-distance part of the integral estimated with the leading OPE 
behaviour, the resulting expression for Bk is 



B 



K 



1 



Mhad 



32tt 2 F 2 



a S (^had.) 

ay log 



a s (rtiad.) 



1229 



z + pv 
Pv 



fa 



1936 
1 

z + pv Pv 



1 



- log z 

Iv 
2 



O 



A^ c a 2 (/Li h ad. 



1 



{z + Pv? 



Pv 



(4.6) 



As seen in Fig. 1, which corresponds to ^had. — 1-4 GeV, F = 85.3 MeV and My — 770 MeV, the 
intersection of the h adro nic curve and the OPE curve happens at a value z ~ 0.39 (i.e. Q ~ 0.88 GeV), 
at which value Eq. ([4.6| ) gives 

Bk^0.38. (4.7) 

The stability of Bk versus z is rather good, and it is shown in Fig. 2 for the same input values as in 
Fig. 1. 

9 The constant L3 gets also an extra contribution from the lowest scalar state, but it is rather small. 
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0.1 0.2 0.3 0.4 0.5 



Figure 1: Plot of the hadronic function W[z] in Eq. \3.\ ) versus z — QV^had f or Mhad. = 1.4 GeV. 
The red curve is the normalized function W[z] corresponding to the minimal hadronic ansatz discussed 
in the text, with a vector meson of mass My = 770 MeV. The green line represents the low-energy 
chiral behaviour and the blue curve the prediction from the OPE. 



Notice that, if we had had a perfect matchin g be tween the minimal hadronic ansatz and the OPE , 



we could have taken the limit z — > oo in Eq. (4J3) and find again the same result as in Eq. (3 



with the sum over hadronic states restricted to the lowest vector state. The fact that, as seen in 
Fig. 1, the matching is not perfect is not surprising; it is due to the restriction of the infinite sum in 



Eq. (3.1) to just the couplings of the lowest vector state. It is quite remarkable that, already at this 
approximation, the quality of the matching is so good. The advocated choice of a £ which minimizes 
the value of Bk at the separation between a low Q 2 region and a high Q 2 region implicitly assumes 
that the leading term of the OPE controls reasonably well the behaviour of the underlying Green's 
function W[z] down to Q values in the region 0.8 GeV < Q < 1 GeV. 



In order to quantify the errors in Eq. (4.7) we proceed as follows: for every choice of /Zhad. one 
finds the corresponding value of z for which the hadronic ansatz and the OPE intersect. This value of 
z is then used in Eq. ( |4.6| ) to obtain Bk- The error is estimated by allowing for a reasonable variation 
of the input parameters: My — 770 ± 30 MeV and 1 GeV< /Xhad. < 3 GeV. The corresponding values 
for Bk are given in Table 1. They correspond to the spread of values: 

0.33 < B K < 0.44. (4.8) 

One way to quantify the systematic error of the minimal hadronic ansatz approximation of this 



calculation would be to include in the analysis e.g. higher order terms in Eq. (3.3). Until we do that 
we suggest as a cautious rule-of-thumb estimate of this uncertainty to round off the spread in Eq. (4.8) 
to an overall 30% of the central value, with the result 

B K = 0.38 ± 0.11 . (4.9) 

Several remarks concerning this result are in order. 

10 In principle it is possible to impr ove on the minimal hadronic ansatz approximation we are adopting here by allowing 
for couplings of higher states in Eq. ( |4.6| ) provided of course that higher order terms in the OPE (or the chiral expansion) 
of the function W[z] are calculated as well. This is under investigation at present. 
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Figure 2: Plot of B K in Eq. tyjty versus z, for fi had . ^ 1.4 GeV, F = 85.3 MeV and M v = 770 MeV. 
Notice the vertical scale in the figure. 



Table 1: Summary of Bk results for different input values of /ihad. and My 



M v 


^had. 


lGeV z 


2GeV z 


3GeV z 


4GeV z 


5GeV z 


6GcV z 


7GcV z 


8GeV z 


9GeV z 


740 MeV 


0.387 


0.413 


0.421 


0.426 


0.429 


0.431 


0.433 


0.434 


0.435 


770 MeV 


0.388 


0.413 


0.422 


0.426 


0.429 


0.432 


0.433 


0.435 


0.436 


800 MeV 


0.329 


0.342 


0.346 


0.347 


0.348 


0.349 


0.350 


0.350 


0.350 



Our result in Eq. (4.9) does not include the error due to next-to-next-to-leading terms in the 
1/N C expansion nor the error due to chiral corrections in the unfactorized contribution, but we 
consider it a rather robust prediction of Bk in the chiral limit and at the next-to-leading order 
in the 1/N C expansion. 

Our calculation shows a crucial qualitative issue, which is the fact that the low-energy hadronic 
contribution below a mass scale 1 GeV < /ihad. < 3 GeV brings down, rather dramatically, the 



value of Bk\il\i in Eq. (4J2) evaluated at that scale. We can already conclude that any calculation 
which " ignores" the details of the low-energy hadronic contribution will give an overestimate of 
B K - 



• The result in Eq. (4.9) is compatible with the current algebra prediction |2q] which, to lowest 
order in chiral perturbation theory, relates the Bk factor to the K + — > vr+Tr^Mecay rate. In fact, 
our calculation of the Bk factor can be viewed as a successful prediction of the K + — > 7r + 7r° 
decay rate! 

• As discussed in ref. |27j the bosonization of the four-quark operator Qas=2 and the bosonization 
of the operator Q2 — Qi which generates AI = 1/2 transitions are related to each other in the 
combined chiral limit and next -to-leading order 1 /N c expansion. Lowering the value of Bk from 



the large-iV c prediction in Eq. ( |l.S| ) to the result in Eq. (4.9) is correlated with an increase of 
the coupling constant in the lowest order effective chiral Lagrangian which generates AI = 1/2 
transitions, and provides a first step towards a quantitative understanding of the dynamical 
origin of the AI = 1/2 rule. 

Finally, we wish to point out that the techniques applied here can be extended as well to include 
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higher order terms in the chiral expansion. It remains to be seen if chiral corrections to the 
unfactorized term in Eq. ( [2.3[ ) could be so large, (of order 100%?) as to increase our result 
in Eq. (4.9) to the values favoured by the lattice QCD determinations [^] as well as by recent 
phenomcnological arguments p9| . 
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